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ABSTRACT 


Expressions for a generalized Modified 
Transitional 3utterworth-Chebyshev (MTBC) filter are 
derived. The characteristics of this filter as applied 
to digital filter design are investigated. It is shown 
that Бу adjusting location and order of the inserted 
zeros, the cut-off slope rate of the filter can be 


traded for maximum attenuation in the stop-band. 


The perfornance of this MTBC filter is compared to 
that o f Butterworth, Chebyshev, transitional 
Butterworth-Chebyshev filters together with those 
suggested by other investigators [1]-(3]. It is shown 
that the stcp-band attenuation can be significantly 
increased without great sacrifice of cut-off slope 


rate. 


Step response of this MTBC filter is also obtained 
and compared with other filters. Various tabulaticns 
as well as graphs of this filter are given for design 
purposes. A computer program is developed for the 
design of this filter. 
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The digital filter is, as defined by Rabiner et al [10], 
" a computational process or algorithm by which a digital 
signal or sequence of numbers (acting as input) 15 
transrormed into a second sequence of numbers termed the 


Mase dig¢tal signal". 


The area of digital filtering can be divided into two 
major subdivisions as Finite Impulse Response (FIR) filters 


and Infinite Impulse Response (IIR) filters. 


During tke development of digital signal processing, the 
interest cr the investigators in IIR and FIR filters varied. 
Before the introduction of the FFT algorithm by Cooley and 
Tukey (1965) IIR filters were much more efficient than FIR 
filters. Stcecham's work [13] on the FFT method of 
performing convolution indicated that implementaticn of 
high-order PER tilters could be nade extremely 
computationally efficient ; thus, comparison between FIR and 
IIR filters are no lcnger strongly biased toward the latter 
[IS "IA Because FIR filters require very high orders to 
produce a Sharp attenuation shape, they are not often used 
for real-time filtering of waveforms. Recently, due to the 
increase in computing capabilities in digital Signal 
processing and the availibity of long charge transfer device 
(CTD) tapped delay lines(TDL), FIR filters are favored over 
πο TES. However, in applications like design of 


digital ccmb filters IIR filters are the unique alternative. 


There are three basic design techniques of IIR digital 
11155 [5]. 
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First method is the direct design, which is, 
appropriately placing poles and zeros to approximate 


required freguency response. 


A seccnd method is to use an optimization procedure to 
place the roles and zeros to match arbitrary frequency 
response specifications. 

Finally the third technique makes use of highly advanced 
art of continuous filter design. This technigue of 
designing digital filters from continuous filters by means 
of mathenatical transformations is the most popular IIR 


digital filter design technique. 


Standard Z-transforn,Bilinear Z-transform and the 
matched Z-transform make possible direct transformation from 
S-domain то Z-domain, preserving essential characteristics 


of analog frequency response. 


Zxistence of frequency transformations reduces the 
problem to design a frequency normalized prototype low-pass 
filter.Then using appropriate frequency transformation, this 
prctotype may be converted into desired  band-pass, 
band-reject cr high-pass filter. Popular prototype filters 
are Butterworth, Chebyshev, elliptic and hybrid transitional 
filters. Frequency transformations for digital filters аге 


discussed in various literatures ( [6] and [7]). 


The problem of designing low-pass prototype filters, 
which possesses better stop-band attenuation and cut-off 
slope characteristics than existing prototypes has always 


attracted the researchers in the signal processing area. 
Budak and Aronhime suggested [1] modification of 


Barimally flat rational functions by introducing a pair of 


finite transmission zeros such that the maximally flat 
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Characteristic is maintained but the cut-off slope can be 
made steeper wa thou t great sacrifice of stop- band 


attenuation. 


Dutto Roy [2] investigated a more general case allowing 
insertion of multiple pairs of transmissicn zeros, either 


Gommerdaent or distinct. 


Introducing multiple pairs of jw-axis zeros in all pole 
Chebyshev transfer functions are investigated by Agarwal and 
sedra |3]. 


The most attractive feature of these finite zero filters 
is that they offer the filter designer a great degree of 
freedom in choosing the location and order of the zeros to 


trade cut-off slope for stop-band attenuation. 


In this thesis, a nodified Transitional 
Butterworth-Chebyshev filter is developed, which is a more 
qeneral case, introducing finite  coincident or distinct 
multiple pairs of transmission zeros in transitional 


Butterworth-Chebyshev filter. 


Trade-orf's between the order of the filter, the order 
of transmissicn zeros, stop-band attenuation and cut-off 
Slope are pointed out. Graphs helpful in the design of such 


Filters are crtained. 


Performances of Butterworth, Chebyshev, Transitional 
Butterworth-Chebyshev filters and the designs suggested in 
references [2] and [3] are compared with those represented 
in this thesis for the orders of three through eleven. A 
computer program is developed to implement the filters 


mentioned abcve. 


In addition, the time-domain response of digital filters 
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is studied. “There are many applications, such as digital 
MTI filters, for which one is interested in the transient 
responses cf filters that are specified in the frequency 
domain. Stef responses of the filters that are discussed in 


this report are plotted and compared. 
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II. DERIVATION OF MODIFIED TRANSITIONAL 
BUTTERWORTH-CHEBYSHEV FILTERS 


А. INTRODUCTION 


The most popular technigue for designing IIR digital 
filters is to digitize an analog filter that satisfies the 
design specifications [5]. There are many techniques for 
designing analog low-pass prototype filters. Among the well 
kncwn analcg IO) Uer classes are the maximally flat 


(Butterwcrth) and equal ripple (Chebyshev) filters. 


Butterworth filters are simple, excellent in the 
pass-band and monotonic in both pass-band and stop-Lband. The 
Chebyshev filters are superior at and near cut-off frequency 


amd at stop-rand. 


The transitional Butterworth-Chebyshev (ТВС) filters 
СРБЕ Ое desirable attributes of these two filters іп а 


single appgrcximation that is given by 


5 1 
ғ = aya CG ey 


where 

k= Weighting factor 

n= Order cf filter 

E (n-k)** order chebyshev polynomial. 
When k=n, ІЕ (3ҹ) {2 is identical with the Butterworth 
function. when К-0, IF(jw) |2 is identical with the Chekyshev 
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function. With a varying value of k, a TBC filter rossesses 
scme characteristics of each. As k approaches to n, a TBC 
filter behaves more like a Butterworth filter, as k 


approaches to zero, it behaves more like a Chebyshev filter. 


ντ ταις Chapter, modification of TBC filters by 
Zmeroducing finite coincident or distinct multiple pairs of 
transmission zeros will be discussed. It will be shown 
that, using a weighting factor k, and the location and order 
of inserted zeros as parameters, attenuation in the 
stop-band may be traded for sharpness of the cut-off 


Characteristics. 


Expressicns £or the cut-off slope and minimum 
attenuaticn in the stop-band are derived in terms of order 
BN he τιίτθας, weighting factor, location, and order of 


inserted zerces. 


Poe err GAM TON OF THE TBC FILTERS WITH COINCIDENT 
TRANSMISSION ZEROS 


Introducing m identical pairs of transmission zeros at 


zu to eg. (2-1), we have 
2 (ως - wt)” 
К бу) | BELLNM Y ise 52 y. (2-2) 
AA AA E 


In order to normalize |F(jw)|?, i.e. to force |F(jw) |? to be 


equal to 1/2 at w-1, the constant K should be 


-—— 4) 
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Then eq. (2-1) becomes 

| Fl š ) Γ G с ыз) sm @ 4) 
u) Ξ- ا‎ = 

) (QR iy^ E Сюе 1)?" “یں‎ QUNM Cw) 


where n>2m, Eecause of the low-pass characteristics of the 


function, and the Chebyshev polynomial C2 (w) is defined by 


Cos” (n Cos w) > wl (2-5) 
C. бо) = 
Cosh” (n Osk w) z ία Т; 71 (2-6) 


¡silo pe at Cut-off frequency 


SubStitutilng eq. (2-5) into eq. (2-8), we obtain 
Am 


|F | MN ο; 
Cu) 7 CY w*) EM ο η uy Z Cost ΚΕΠ δος! ιο] 


L 63) 


f(w)= (we ut)" (2-8) 
h (10) = ае Cos” ΚΟΠ Cos tw | (2-3) 


Using these values, eq. (2-7) becomes 


2 + (w) : 
| FG) | “yr MEA πίω). D 
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Taking the derivative of eq. (2-10) with respect to w, we 
get 


ο (iio [C3 G) - 662 k (J 
КӨ K 2 |FG»)| ЕСЭ Сах ο (4-9) 


At the cut-off frequency, wz1, we have 
Flo) = (wast 4) (ег) 
αμ) La T (2-13) 


№ (о) = 4 (2-14) 
h'(w)= 2k + 2 (n-2)? (2.15) 


IFGw)| = 1/11 (2-16) 


Surstituting these values into eq. (2-11) and simplifiying, 


we obtain 


ki (n-k)* 
FG) | = | ya Cor Ы "Жұл ер 


For k=n, the result agrees with the cut-off slope of MB 
cein, Which is derived in [2]. For k=0, the result 
agrees with the cut-off slope of MC function as derived in 


L3]- 


SES band characteristics 


-— OA C πο “as ma s ны = Sa á- e M e å сс å эшш πα x m= сс 


In general, stop-band characteristics of finite zero 
ЕЛІсСерве Mill re of the form shown in FIG.1*®. 


* Because of the inherent limitations in the plotting 
әншроцріпе utilized to provide the graphs in this thesis, it 
was necessary to add supplementary axes to show proper 


scaling for scme of the graphs. 
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The zeros introduced will cause a peak in the stop-band, at 
a frequency Wip Wo . The minumum attenuation in the stop-band 
may be defined as 





ο ρε = 20 Log To (2-13) 
where 
Te = Р] (2-11) 
In the stop-tand eg. (2-4) becomes 
2 τ. t u) aver 
> | = D. (n-)- 2 (г) 
(м P) “+ Cw Dd i De E ust 
Let 
С) = Cae. vi Е ! (2-21) 
2(n-k)- 2n. 
h ζω) = 2 Wr (2-22) 


Taking derivatives of these values and substituting them 


ЮП О ец. (2-1) results in 


Fi = 
| (up) | 4 [К (фә) | Κο e ,40092 ui D 


Combining 29. (2-19) and eq. (2-23), we obtain 


We = . ἘΝ (2-24) 


n- Am 


— n E 
A S $29 m-k-1) n m) 
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And minugum stop-band attenuation will be given by 


Xue = 100 Log EX +20 m 109 (Wo) 


+20 m Log nim WA - 6 (n-m-h-1) (2-26) 


Wi 


Plots of stop-band attenuation and cut-off slope of 
EUM filters with two coincident transmission zeros, for 


oarders 3 thrcugh 11 are given in figures 2 and 3. 


ο οὉτεκιτοουνοντνς FILTERS WITH DISTINCT TRANSMISSION 
ZEROS 


th 


Consider the a order TBC function with a pairs of 


nee distinct zeros attjw , where i=1,2,...,m and vw; 21. 


Phe @#adgnitude squared function of TBC function with 


distinct transmission zeros will be given by 


ғо) |“ = Л Goa (ἢ 
u) — N —əə - n 27. 
р Toot HOSE ë Te tay oe Cy (a) 


Let 


ІМ 


alw) = Т «Құлан 
πίω) - ; Τζων. ων) 


— these values into eq. (2-27) resuits in 


2 
ғы) | + | Aw) — ο ο) (2-23 
6 πο 


22 


E 














EL 


ΠΑ. Ju DL“ 


اك 


2 5 


stop - band attenvation CU) um 








. Ne eee O ο πο ы 
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2 / 
IE 
Ë / 
- П m y ae EX ^ 
4.2 "s 2.9 2 3.0 7.5 4.0 


Figure 2 - STOP-BAND ATTENUATION OF MTBC FILTER WITH TWO 
COINCIDENT ZEROS 
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ους Band ald cut-off characteristics 


a > сш о шшш сы сш ع‎ qam κο Ñ m t s = u 


The frequency wp at which the stop-band peak will 
occur, is given by 


[FG |. - 0 (223) 


م لہا <= لا 


Using eq. (2-6) and taking derivative of eq. (2-28) with 


Meepect tc w, we obtain 


| Elsa | ' š аш) ο τα k )-2 n lan πω) ТІ! 05) ] (2:30) 


G@em@bining eq. (2-29) and eg. (2-30) the equation to he 


Sclved for Wo may be found to be 


-- 221 
2 We $ "i πο |... O ( ) 
(=) (ωιτ-ωρ) 


which agrees with eq.(2-24) when w;'s are identical. Sang 
real soluticns of eq.(2-31), stop-band attenuation peaks 


may be rcund to 5e 
X, = 20 log Tei (232) 


where 


To; = | F Gg) (2.23) 


and minimum stop-band atttenuation is given by 


Pd c mar TT (2.24) 
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NEMUS eq. (2-5) and differentiating eg. (2-27) with respect 
to w, the cut-off slope becomes 


m 
FG) | E = R + (n-k)* + el P we (2.35) 


which agrees with eq. (2-17) when the w, 's are identical. 


Plots of cut-off slope and stop-band attenuation of 
ВС filters with two distinct transmission zeros, for 


GEES 3 thrcugh 11 are given in figures 3-17. 
D. SUMMARY 


The Transitional Butterworth-Chebyshev filter combine 
the best features of the Butterworth and Chebyshev filters. 
A Modified Transitional Butterworth-Chebyshev filter, 
obtained by introducing finite transmission zeros to a 
Transitional Butterworth-Chebyshev filter, possesses cut-off 
slope and stop-band attenuation, which are dependent on the 
modification parameters wo and m. The closer the ws is to 
unity, the steeper the cut-off slope. However, this 
improvement in the cut-off region results in degradation of 
stop-band attenuation. Thus, using w, and m as parameters, 
an advantageous trade between attenuation in stop-band and 
Sharpness of tne cut-off characteristics can be made. 
Graphs are given to serve as guides in trading cut-off slope 


for stop-tand attenuation. 
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A. INTRODUCTION 


ο μα τες 2, derivation ot MTBC filter is given. It is 
Shown that, stop-band attenuation and cut-off slope rate of 
this filter depend on modification parameters m, k, and we. 
So the filter designer have the flexibility to trade cut-off 
slope for stop-band attenuation Dy changing these 


parameters, without changing the order of the filter. 


In this chapter, performance of MTBC filter is compared 
MEN В, ME, C, MC, and TBC filters. 


Formulas for stop-band attenuations and cut-off slopes 
of all these filters are given in Table I. 


БОКУ ИНС MONCTION V.S. B FUNCTION 


οἱ off slope of MTBC function is given by 


E m k(k+-2n) I 
nac = E Ya" d ae m € 1) 
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ου ο ratio of cut-off slopes of these two 


functions may be found as 


ED JZ m πώ ου 
ἵμτας /6 = тє - 1 + nO) pO ur (3-2) 
T seca V.S. Ww, for n-3 through 11 and m=1,2 are given in 


quese 18 and 19. Stop-band attenuation of MTBC may be 
written as 


Aura = gr 6(n-m-k-1) + 20 m Log A EX = e 9) 


Then the difference between stop-band attenuations is given 
by 


аа XK тес - оса = 6 (n-m-k-t) +20 m log IA. E 


(4) 
Plots of άμτος/α Ves ο τον through 11 and m=1,2 


are given in FIG 20 and 21. 


Figures 19-21 indicate that, MTBC filter can be made 10 
times steeper than B filter still having 40 dB more 


attenuation at tne stop-band. 
CONMNERBBC U.S. WB FUNCTION 


Melina Table I, cut- off slope ratios of MTBC and MB 


filters may be found as 


СЕ | | (n-k) eke | 
nrac/me = 2m +n (we-4 (5-6) 


БІНЕ Of r V.S. W for n=3 through 11 and 1-1,2 


mtac Ju 8 
given in figures 22 and 23. 


Tne difference between the stop-band attenuations of  MTBC 
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and MB is given by 


тес (ав E 6(n-k-4) 


Baq@ation (3-7) indicates chat; stop-band attenuation 
difference between these filters doesn't depend On Wa. For 
n=10, m=1, k=1, with 48 dB more attenuation, cut-off slope 
ο ο жау τε changed from 2.5 to 8.0 by changing Wọ fron 
HIJOS to 2.4. 


E TEES FUNCTION .¥.5. C FUNCTION 


Mie “ratio Of Cut-off slope of MTBC function to cut-off 


slcpe of C function is given by 


Im k.(k+1- 24] 


Γ ο. ow munt i Lp 
мтвс/с V asl) Ai 


Plot of T urge /c V.S. w, for n=3 through 11 and m=1,2 are 
given in Figures 24 and 25. The difference between 
ED band attenuations of MTBC and C functions is given by 
E 
n-2m  Wo”-1 
тве fe = —6(mxk) +20 m Log [tm Wo” | 
Plots of 3urec/c Wes We for n=3 througn 1! and m=1,2 are 
given in Figures 26 and 27. The Chebyshev filter is known 
to provide much steeper cut-off slope than the corresponding 
B and TSC filters. Figure 25 shows that the cut-off slope 
of the MTBC filter can be made 1.4 times steeper than that 


of the Chebyshev filter, with n=5, m=2, and wo=1.06. 


πε FUNCTION V.S. EC FUNCTION 
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me ratio of cut-off slope of MTBC filter to cut-off 
slope of MC filter is given by 


onde EC e klk 2n) (ω 1) EN 


Фар Қ Αα ζωσ-{ ) 


Plets of Г атес/мс V.S. W for n=3 through 11 and п-1,2 аге 
given in Figures 28 and 29. 
The difference between stop-band attenuations of MTBC and MC 


functions is given by 


Сыз ας --6k (2-11) 


For a given k, c is allways constant. 


Û rec /u 


fo meee Y.5. “EBC FUNCTIONS 


Mae parao. of the cut-off slopes of MTBC and τες 


Bumictions is given by 


Im 
pr A een. — 3-11 
τες /τες + TEN [nta (841-4431 ( | 


Pots от Гитас/тас "<S. We for n-3 through 11 and m-1,2 are 


given in Figures 30 and 31. 


The difference between the stop-band attenuaticns of 


these two functions is given by 


үн = — 6m + 20 m Log | s | N (2-13) 


procs of dutec/tec VeS. Wo for n=3 through 11 and m=1,2 


аге Чеп in Bigures 32 and 33. 
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б: SUMMARY 


Using the location and the order of the inserted zeros 
and the weigkting factor as parameters, the characteristic 
curve of the MTBC filter can be made steeper than that of 
the conventional all-pole filters without greatly 
sacrificing either stop-band attenuation or flatness in the 
pass band. A Modified Chebyshev filter provides slightly 
better performance than the MTBC filter, at the expense of 
substantial degradation of the pass-band flatness. Graphs 
are given to help in the comparison and in determining the 


numerical advantages gained by increased complexity. 
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А. INTROLUCTION 


in the previous chapters, expressions for MTBC filter 
are derived and the performance of this filter is compared 


Were various filters. 


In this chapter, A computer program is developed to 
¡aplement MTEC filter. ten ishi pointed out that all five 
filters, which are used to ccmpare with the MTBC filter, are 
the special cases of the MTBC filter. Thus the program 
developed in this chapter for the implementation of the MTBC 
filter may also be used to implement any one of these rive 


БІНЕ ӨТ. 


À sample problem is worked out to illustrate the use of 
the plots presented in earlier chapters, the use of the 
computer  prcgram, and to point out the flexibility offered 
OT he το filter to the filter designer. 


B. COMPUTER FROGRAN 


In the most general form the transfer function of MTBC 
ENIM -cuJven by eq. (2-27), which is repeated here for 


T (ωύ- ω”)” 
T (wi? Y. ‚Ta 4-4) o** бқ бы) 


e 


convenience. 


اد 





a 


Fo] = 


where : 


. 4h 


= = Location of i inserted zero 


Ka Order of inserted zeros 


Bm = Order of filter 


k = Weighting factor of Transitional Butterworth 
Chebyshev filters 


When k=n, and m=0, |Е(34)|2 іс identical to the 
Butterworth function. With k=n and m#0, |F(jw) {2 is 
identical to the modified Butterworth function which is 
discussed by BUDAK and ROY [1] - [2]. When k=0 and n=0, 
iF (jw) ]* is identical to the Chebyshev function. With k=0 
and m#0 Modified Chebyshev function may be implemented which 
is discussed by AGARVAL and SEDRA [3]. 


When m=0 and τεῦ, ІР (9ҹ) 12 is identical to the 
Transitional Eutterworth-Chebyshev filter. Finally allowing 
both k and m το vary, Modified Transitional 


Butterworth-Chebyshev filters may be implemented. 


ШЕ ргсагал consists of two main parts. Іп the first 
БЕС СИ analog filter, which is given by eq. (2-27), is 


izplemented and the frequency response is plotted. 


In the second part, the analog transfer function F (jw) 
is first predistorted then transformed into z-domain by 
algebraic substitution method (using Bilinear 
Z-teaiaseermation) to obtain H(z). Digital transfer function 


{ H(z) ] is then factored into second order cascaded stages. 
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“ΤΥ, freguency respoase curves ( q|F(jw)|. v.S. wy 20 


loq |F(3j<q) | v.s. wand Fe ws. Wo.) are drawn. 


Bilinear Transformation is prefer 
available algebraic substituticn a (246% Impulse 
invariant, Watched z-transformation) in obtaining E(z), 


ily ісе the following reasoas [4] : 


(1) it nas the btoperty that realizable stable 


continuous systems are mapped to reaiizabie stable digital 
š 


(2) Wideband sharp cut-off continuous filters can ре 
Берре to »sbedeoand sharp cut-ofz digital filters without the 
aliasing in the frecuency response. 

(3) After Bilinear Transformation, the relation between 
tue analog amd digital zrequencies is given by 


A 
By choosing cut-off frequencies approaching to fs/2 where És 
Stands ior sampling frequency, extremely sharp cut-off 
u 


ers gam Be obtained. Ms an examPie plot of c 
V 


Of bilinearly transformed Butterworth filter 
3i 


4 


frequency is given in FIG. 


BATA CARDS 


Q 
LU 
t 
O 
C 
ні 
j 
tri 
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Tuc Wadia cards required to use the program are given 


below. 


ӨШЕй 77 : Values of a,a,k in 3I2 foréat 
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νη: Location of inserted zeros, и; where i=1,2,...,m 
in 8F10.5 format 


Card 3 : Initial and final values of the frequency to be 


used for the frequency response plot in 2F10.5 format 


Card 4 : Number of solutions required in I3 format 


DT REVUIRED SUBROUTINES/FUNCTIONS 


In addition to the built-in subroutines, the following 


IBM source library subroutines are used. 


Ш- POL ET 


4. ΕΡΕ 


9. Powe 


q. PHEY 


5. PADDM 


ЕБ; DESIGN £XAMPLE 


Suppose we want to design a digital filter with a cut 
Sie jO COE 15 ama winumum attenuation in the stop band of 
60 Db. 


Pass lle solutions for various types of filters are 


given in table II. 
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SLOPE \ATTENUATION 


reater than 


MB Lo 





| 
x 
| 





|9045 
15.05 (ДЕР 
rea ter Shan 
6 5 
60 20-8 | 
| 
20.5 60 78 
| Greater 748/7 
ЕЙ ТРЕ A d e 





PNIS C 


EL 


=— = —— ο P © 


3 


1 


ЛА 


7, /f 


69 


| نے سے د‎ — |——— -ϕ 








Enc motlowmsng input data specifies a МТВС filter of 
order 7 with m=2, w,=1.36, k=1, sampling period of 1 sec., 
100 soiution points, and frequency response plot from 0 to 3 
BZ. 


7181 31191 11 [12113 [14115116117 |18 1 19 1 20121 1221231241251 261 271281291 30 


Frequency responses of three of the possible solutions 


Or Ehe design example are given in Fig.35. 
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V. TIME DOMAIN RESPONSE OF DIGITAL FILTERS 


— mmm сс = == == = m с mm Qan EEE mb m =s s Qs sm ae ج‎ самас چ‎ COE ms mtb s s m ιν 


А. INTRODUCTION 


In previous chapters we have investigated the problem of 
desining digital filters in the frequency domain, i.e. to 
meet given frequency domain specifications. A filter 
designer should always consider the transient response 
characteristics of its filter. There are many applications, 
Such as digital MTI filter, for which one is interested in 
the transient responses of filters that are specified in the 
frequency-domain. Tine-domain and frequency- domain 
characteristics of a filter will work against each other. 
Filters close to the ideal frequency characteristic can be 
designed. Filters whose time characteristic is close to the 


ideal can also be designed, but filters close to both cannot 


ie. 


In this chapter time~domain response of digital filters 
will be discussed and it will be shown that the location of 
the transfer functions poles has a profound effect on the 


transient response of the filter. 
B. TRANSFER FUNCTIONS! POLES AND TRANSIENT RESPONSE 


EN а transfer function of a digital filter іп 


factored form 
ЩЩ (ъз). tt) | 
H (2) = ХО (ар) (2-р): (Pm) (5 1) 
7 


2 








where Z,'s and pj 5 are the zeros and poles, respectively, 
of the filter. In general, the input signal's Z-transform 
So: the form 


б. N (z) i 
a (2-4,) (8-42) G=qe) (ra) 


where qi's are the poles of the input function. The 


response of the system may be obtained from the transfer 


function relationship 


Ya) = (а). (а) = EI em) |] NG) | 
G-A) (а-а) (2 -q,)-- (2-4) 
(553) 


PO” Slimplicity, we assume that all of the poles of the Ү(2) 
BEN distinct, and making partial fraction expansion of Y(Z), 


we obtain 


p, z Rez C, 2 Gace 
2) = Ro + E WE. ens mcr as — Ñ ... А 
Ү( 2-4. Po + is + ἘΠ +: + сз 


(5-4) 


The response of the system may be decomposed into [5] two 





parts called the input signal node (IS4M) and the system mode 
(SA) as 








к. Riz R, + Era 
CAE бе ne 5-6 


Combining eg. (5-4) with eq. (5-5) and eq. (5-6) 


Y (a) Eds омы >) + Ya (29) (5-7) 


Eq. (5-7) indicates that the response of any linear system 


to any input will contain modes generated by the input 
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B m = 
Т ы .-- 


4 
o 





Signals roles апа the system transfer function poles. Тһе 
time domain response of the filter may then be found by 
taking the inverse Z-transform of eq. (5-7) 


yk) = Re dk + Yisulk) + Jsu Ck) (5-2) 


where 


Waq (k) = k, (409 ka CA "ELEC E @) 
Ysu Ck) = б, (ру πο Сы) E Cg.) (5-10) 


my@attons (5-9) and (5-10) indicate that a fast responding 
System is one in which all of the system transfer function 
HESSE pare sufficiently smaller than unity in magnitude, 
in order that the system mode will decay to zero rapidly. 
On the other hand, a slowly responding system is one in 
which the system mode decays to zero very slowly (i.e. at 


least one of the p is close to unity in magnitude). 


Those poles closest to the unit circle will be called 
the system's dominant poles because they tend to dominate 


the characteristic of the resultant transient response. 


The dominant poles' response property is, as stated by 
Cadzow [ 14] "The response tine of a linear discrete system 
is directly dependent on the locations of the systen 
transfer  function's dominant poles. Depending on the 
particular response-time requirement FOE a given 
application, we then have to correspondingly locate the 
dcminant poles of the transfer function. A fast responding 
system necessitates dominant poles of magnitude much less 


than one ", 


The poles of the MTBC filter for the orders 3-11 and for 


various values of m and w, are given in Table III. 
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Although this partial fraction expansion method helps us 
to understand the importance of the poles of the systems 
transfer function in transient response analysis of digital 
filters, evaluation of residues of corresponding poles of 


partial fraction expansion is not a trivial problem. 


we believe that the so called ' transfer matrix " 
method, which we are about to discuss, is more suitable for 


digital computer simulation. 


Given a digital filter weighting sequence h(n) and the 
input sequence x(n), the response of the digital filter may 


be obtained by convolution summation | 


y Un) = 2 h(n-k) x Ch.) = y hlk) x (n-k.) (=u) 


Eq. (5-11) may be written in matrix form as 






y (0) hlo) ο ο. ΠΠ X93, 

yla) kG) Мо) 5. 

: = | : ; . : (5-1) 
[IA] [μα kan κίοὴ-:. А] 


Or 


У = G. xX б-а) 


where G is the systems transfer matrix, which is defined as 


h Lo) о Он Мек» С 
κα) πο ον ο ο 


(5-14) 


IO 
I! 


Мк) hbr-1) hla-2) --- hlo) 


To find the systems response using eq. (5-13), systems 


transfer matrix must be available. To obtain the weigting 


Το 





Sequence (impulse response) of the system to form the system 
transfer matrix in terms of the transfer function's 
coefficients, the difference equation of the system is 
solved ror the impulse input, i.e., 


Na) а. _ жу bi y (n-i) O (5-15) 


е 
= | 


Maca: ава b- are the numerator and denominator 
e ficients, respectively, of the transfer function of the 
digital filter. 


À ccnputer progran ( FORTRAN ) is developed to 
investigate the transient response of the digital filters, 
using equaticns (5-13) and (5-15). Program listing is given 


in Appendix B. 


Step responses of MTBC filter for various values of n,. 
Ku sad ES are given in figures 36-39. Figure 36 and 37 
indicate that increasing the values of nm and wg also 
increases the overshoot and settling time, but doesnt have 
any significant effect on the rise time. The rise time tends 


to increase with increasing order of the filter. 
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Transient responses оё digital filters depend on the 
ВЕЕР оп ов the poles oi its transfer function. & fast 
responding filter has poles of magnitude much less than one. 
Modification of dll pole filters increase settling time, 
decrease peak overshoot, and doesnt significantly affect the 
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rise time of the filter's step respons 
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